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' Spectral functions relevant in the context of quantum field theory under the influ- 

ence of spherically symmetric external conditions are analysed. Examples comprise heat- 
kernels, determinants and spectral sums needed for the analysis of Casimir energies. 
First, we summarize that a convenient way of handling them is to use the associated zeta 
function. A way to determine all its needed properties is derived. Using the connection 
with the mentioned spectral functions, we provide: i.) a method for the calculation of 
heat-kernel coefficients of Laplace-like operators on Riemannian manifolds with smooth 
boundaries and ii.) an analysis of vacuum energies in the presence of spherically sym- 
' metric boundaries and external background potentials. 

> : 

1 Introduction 



There is a notorious appearance of spectral functions in many branches of mathematics and 
physics. These functions are associated with suitable sequences of numbers {^fc}fc g ]N' w hich 
for most applications are eigenvalues of certain interesting, in most applications Laplace-like, 
operators. A rich source of problems where spectral functions are encountered is quantum 
field theory under the influence of " external conditions" . External means that the condition 
is assumed to be known (as a function of space and time) and only appears in the equation 
of motion of other fields, which are to be quantized under these conditions. Given that the 
focus of interest is on the influence of the external conditions, these other fields are often 
assumed to be non-selfinter acting. It is in this setting that the present work is to be viewed. 

Under the described circumstances the action of the considered theory will be quadratic 
in the quantized field. Using a path integral formulation to describe the underlying quantized 
theory, one encounters Gaussian functional integrals which lead to functional determinants 
formally defined as Ilfc^i ^k- Making sense out of this kind of expressions is one of the basic 
themes of quantum field theory. The present contribution deals with scalar fields and will 
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develop and apply techniques to analyse determinants arising when the external conditions 
are described by boundary conditions (Casimir effect) or when external scalar background 
fields are present. We restrict to scalar fields because the application of the provided ideas 
to spinor fields or electromagnetic fields is then immediate (see f.e. [Q, ||). 

In some more detail the situation considered will be described by the action (we assume 
an Euclidean formulation) 



S[^} = ~\ I dx$(x)(n E -V(x))<f>(x), (1.1) 

2 JM 



for a scalar field <I> in the background potential V{x). Here, M is a D-dimensional Riemannian 
manifold and dx its volume element. 
The corresponding field equation is 

{U E -V{x))*{x) = 0. (1.2) 

If boundaries are present, the equations of motion are supplemented by boundary condi- 
tions to be specified later. 

Physical properties of the systems are conveniently described by means of the path-integral 
functionals (an infinite normalization constant is neglected) 

Z[V] = f D$e- s W, (1.3) 



where the functionals are taken over all fields satisfying, if applicable, the boundary conditions 
imposed. 

Under the circumstances described, the effective action can (at least formally) easily be 
computed to be (assume that there are no zero modes; if these are present they have to be 
omitted because otherwise the determinant is trivially zero), 

F[V} = - lnZ[V] = - In det [(-D E + V(x))/tf 

Here, fi is an arbitrary parameter with dimension of a mass, to make the argument of the 
logarithm dimensionless. The operators involved are thus Laplace-like operators, extensively 
dealt with afterwards. Let us write them in the unified form also used later, namely 

P = -g pu V p Vu - E, (1.4) 

where g pu is the Riemannian metric of the manifold M, V is a connection and E an endo- 
morphism defined on Ai. We are thus confronted with the task of calculating expressions of 
the type 

r[y] = iln[det(P/ M 2 )]. (1.5) 
Clearly, expression ( |1.5[ ) is not defined because the eigenvalues A n of P, 

P(f>n = K<t>n, (1-6) 

grow without bound for n — > oo. Of course, there are various possible regularization proce- 
dures; let us mention only Pauli-Villars, dimensional regularization and zeta function regu- 
larization. Here, we will use zeta function regularization 0, 13 because it is mathematically 
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appealing as well as (probably) the most convenient one in the context of our work. The 
basic idea is to generalize the identity, valid for a (iV x iV)-matrix P, 



N ^ N , 

lndet ^ = E lnA ™ = E X n S \s=o = -^Cp(s) 

rt=l n=l 



with the zeta function 



N 

Cp(*) = E A n S , 
n=l 



to the differential operator P appearing in (|1.6| ) by 

lndet P = -Cp(0), (1.7) 

with 

oo 

CpW = E A n S - (1-8) 

n=l 

That this definition is in fact sensible is a result of deep mathematical theorems on the 
analytical structure of Cp( s ) (explained in the following). 

First of all, due to a classical theorem of Weyl ||, which says that for a second order 
elliptic differential operator the eigenvalues behave asymptotically for n —* oo as 



D/2 2 D -^DT(D/2) 
n ~ vo\{M) 



n. 



the representation ( |1.8|) of Cp( s ) is valid for 5Rs > D/2. In order to use definition (|1.7| ), the 
question arises of how to analytically continue Cp( s ) to the left and to determine its analytic 
structure. This is very elegantly done by using an integral representation of the T-function 
to write (see f.e. @), still for 3f?s > D/2, 

1 poo 

CpOO = / (1-9) 
r s Jo 



with the heat-kernel 



Jf(t) = E e 



n=l 



For t oo the integral is well behaved due to the exponential damping coming from K{t) 
(for simplicity we assume a positive definite operator P). Possible residues only arise from 
the t — > behaviour of the integrand, and we need information on K{t) for t — > 0. This is 
given by the heat-kernel expansion 0, ||, || 

oo 

tf(t)~ E az(^ C/2 , (1-10) 

£=0,1/2,1,... 
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with the heat-kernel coefficients ai(P) depending, of course, explicitly on the operator P 
together with the boundary conditions chosen. If the manifold has no boundary, the coef- 
ficients with half-integer index vanish. Splitting the integral for example into Jq 1 dt + f£° dt 
the following connection is obtained 



Res (Cp(s)T(s))\ s=d/2 ^ = ai (P), (1.11) 

or, showing the information contained more clearly, for z = D/2, (D — l)/2, 1/2, — (2n + 
l)/2 5 neM , 

Res Cp(z) = ^'IZ^ ^ (1-12) 

and for q £ INo, 

Cp(-q) = {-l)U\a D/2+q {P). (1.13) 

This clearly shows that, for manifolds without boundaries, the poles are located at z = 
D/2, D/2 - 1,...,1 for D even, and z = D/2, D/2 - 1, 1/2, -(2n + l)/2, n e JN , for D 
odd. In addition, for D odd and q € INo, one gets (p(— q) = 0. For manifolds with boundary, 
additional possible poles appear and here one finds poles at D/2, (D — l)/2, 1/2, — (2n + 
l)/2, n £ INo- In all cases, Cp( s ) 1S an analytical function in a neighbourhood of s = such 
that eq. ([O]) is well defined. This definition was first used by the mathematicians Ray and 
Singer |ll] , when trying to give a definition of the Reidemeister-Franz torsion [^] . 



Let us stress that the expansion ( 1.10 ) of the heat-kernel and the connection ( 1.11 ) with 
the zeta function strongly depend on the assumptions made, and strictly they hold only if we 
are dealing with a second order elliptic differential operator on a smooth compact Riemannian 
manifold with a smooth boundary and local elliptic boundary conditions 0. For example, 
in the context of admissible pseudo differential operators and global spectral boundary con- 
ditions, different powers of t can be involved in ( |1.10| ) and, in addition, (lni)-terms might 



appear. This was noted starting with [13| and later, for example, in [14, [15|, [L6|], where the 
asymptotic expansions for these cases can be found. That indeed generically all possible 
terms are present (what means generically can be made precise) has been shown recently 



[17]. As a consequence, the simple poles of the zeta function may well be located at other 
points depending on the power of t appearing in the asymptotic small-t expansion of K(t), 
and, if there are lnt-terms, double poles may be present as well. Also in this more general 
context, analogous relations to equation ( |1 . 1 1| ) may be stated [18], but here we will make no 
use of these and, therefore, we won't bother to state more details. 

A case of particular interest is the manifold M. = 5 1 x M s with the operator 

P = -^ + P s , (1.14) 

and P s Laplace-like. Imposing periodic boundary conditions in the r-variable this is finite 
temperature quantum field theory for a scalar field and the perimeter (5 of the circle plays the 
role of the inverse temperature. Under the assumption that the potential V and the boundary 
conditions are static, P s is a purely spatial operator and depends only on coordinates on M s . 
In this context, it is known that the energy of the system is by definition 

B = -£ M = -5Sj5<W< >' 
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and it seems natural to use 



E vac = lim E = -FP Cp„(-1/2) 

p— >oo Z 



1 



2V4?r 



r,2 



;i.i5) 



with the scale /2 = (fie/2) as the definition of the vacuum energy. Here, FP means finite 
part for cases where Cp s (— 1/2) has a pole. Due to the presence of the arbitrary scale it 
is seen that the vacuum energy is ambiguous, which generally causes problems to achieve a 
physically sensible answer. However, this is the way the Casimir energy is usually defined 
(see for example pT 



a 24, 1|, 25, ||, 27, |28, M, M]), the idea for the derivation 
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going back already to Gibbons |31[|. The ambiguity may be discussed within renormalization 



group equations found by demanding |3^, 33 1 



dfj, 



0. 



Eq. ( 1.15| ) then clearly shows that the total energy of the system must contain all terms 
present in dD^iPs) i n order to define the needed counterterms and running coupling con- 
stants. These terms describe the energy of the external fields or are the energy needed to 
get a model for the boundary conditions. If cld/2{Ps) 7^ 0, the Casimir energy is determined 
only up to terms proportional to CLD/iiPs) an d this finite ambiguity can (in principle) only 
be eliminated by experiments [28]. A possible way to fix the ambiguity for massive fields will 
be discussed in Section f|. If, on the contrary, ajj/2(P s ) = 0, eq. (|1.15|) gives a unique answer 
for the energy. 

Eq. (1.15) is also the definition one would naively use. To see this, write the Hamilton 
operator formally as 

H = Y J E k (N k + \ 



with iVfc the number operator, to obtain for the vacuum energy 



En 



1 



:<0|tf|0 >= -Y, E k- 



;i.i6) 



The regularization 



E,, 



\s=0 



<p.(s- 1/2)1 



s=0 



i, 2s ;; 2 

l -FP Cp.(-1/2) + \ Q + ln/r j R(<s ,/..(- 1/2) 



k 



1 



FPCp,(-1/2) 



(1.17) 



is clearly equivalent to eq. ( |1.15| ) , the only difference being that renormalization now involves 
infinities. 

Let us stress that, although we have focused on non-selfinteracting scalar fields, the cal- 
culation of determinants like in eq. ( |1.5| ) is of great relevance in several areas of modern 
theoretical physics. In many situations the potential is provided by classical solutions to 
nonlinear field equations. Quantizing a theory about these solutions leads to the same kind 



5 



of determinants as discussed above p4j| . The classical solutions involved may be monopoles 
mUl, sphalerons || or electroweak Skyrmions ||, ||, g(| ||, ||, ||, || g|, g|]. The 



determinant in these external fields enters semiclassical transition rates as well as the nu- 



cleation of bubbles or droplets [47, 48]. In general, the classical fields are inhomogeneous 
configurations and, as a rule, the effective potential approximation to the effective action, 
where quantum fluctuations are integrated out about a constant classical field, is not expected 
to be adequate. The derivative expansion [49| improves on this by accounting for spatially 
varying background fields. Being a perturbative approximation it has however its own limita- 
tions. Having in mind that even the classical solutions are often known only numerically, it is 
clear that it is desirable to have a numerical procedure to determine the quantum corrections. 
Some contributions in this direction are |5(| [5^, |53| . 

To summarize, as we have briefly described, the most relevant spectral functions in quan- 
tum field theory under external conditions are determinants, eq. (L5), respectively the spec- 
tral sum, eq. ( 1.16|) , and the heat-kernel, eq. fll.lCD , describing the ambiguity in Casimir and 
ground state energy calculations. It is the aim of the present work, to provide and apply 
techniques for the analysis of all these spectral functions. As we have seen in the general 
introduction above, it is the zeta function associated with the spectrum which can be used 
as the organizing quantity for all the calculations to be done. 

The starting point of our investigation was the aim to develop a machinery adequate 
to deal with spherically symmetric external conditions, the motivation being that many of 
the classical solutions have this property. In all spherically symmetric problems, the total 
angular momentum is a conserved quantity and, as a result, eigenvalues will be labelled 
by this angular momentum and an additional main quantum number. Compared to one- 



dimensional problems, where many calculations can be done exactly (see f.e. |p4j, 55, 56, |5 



one encounters here the common technical complication of angular momentum sums. Any 
spherically symmetric example can help to understand the difficulties arising therefrom. 

Intuitively, to impose boundary conditions on a spherical shell seems simpler than dealing 
with an arbitrary spherically symmetric external field, because eigenfunctions are known ex- 
plicitly. However, let us stress that analytical expressions for the eigenvalues are not available 
which complicates the analysis considerably. Based on the knowledge of the eigenfunctions, 
an important achievement of our consideration is that various properties of spectral functions 
associated with a spectrum which is not known explicitly can be determined. We explain 
this procedure in detail for the spectrum of the Laplace operator on the three dimensional 
ball || . 

Having a very good (if not complete) knowledge of the zeta function at hand, and us- 
ing the connections explained above, the obtained results are applied to the calculation of 
the associated heat-kernel coefficients and Casimir energies. The applications are grouped 
according to their complexity: i.) the heat-kernel coefficients, which are "local" quantities 
and whose determination is purely analytical, ii.) Casimir energies, which are non-local and 
where additional numerical work is needed. 

For the mentioned reason we start, in Section 3, with the consideration of heat-kernel 
coefficients. An extremely important aspect of this calculation is that it is not just a spe- 
cial case calculation, but very rich information about heat-kernel coefficients for Laplace-like 
operators on arbitrary compact Riemannian manifolds evolves. Supplemented by other (al- 
ready known) techniques involving conformal variations and the application of index theorems 
U, it allowed for the determination of new coefficients for all types of boundary conditions 



[59, 61, g^, |63| (including Dirichlet, Robin, mixed, oblique and spectral boundary condi- 
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tions). 

In Section |] we consider Casimir energies in the presence of spherically symmetric bound- 
aries. Electromagnetic field fluctuations in a spherical shell were considered already a long 
time ago in the context of a simple model for an electron 64, 65, 66, 67], where the zero-point 
energy of the electromagnetic field was supposed to stabilize the classical electron. However, 
as was first shown by Boyer [65], a repulsive force is the result and the simple model fails. 

Our focus here will be on the influence of a mass m of the quantum field, and we will 
consider a scalar field with Dirichlet boundary conditions. Based on the only previously 
existing complete calculation for planar boundaries at a distance R [68], the general believe 
is that for mR 3> 1 the Casimir energy is exponentially small and thus of very short range. 
As we will see and explain, this is due to the planar boundaries and does not hold if the 
boundaries are curved. It might even happen that the Casimir energy changes sign as a 
function of the mass. This surprising result justifies to stand the additional complication due 
to the non-vanishing mass. Although we restrict for simplicity to scalar fields with Dirichlet 
boundary conditions, higher spin fields with various boundary conditions as f.e. MIT bag- 
boundary conditions can be treated as well along these lines (see f.e. [69, 70, |7"lj|). 

Afterwards, we analyze the influence of external background fields [72, |73|]. Although, in 
general, 



not even the eigenfunctions are explicitly given, this knowledge is replaced by the 
information available from scattering theory. So we will express the ground state energy by 
the Jost function, which is known (at least) numerically by solving the Lippmann-Schwinger 
equation. Applying the formalism developed in Section |2] for the treatment of the angular 
momentum sums, quite explicit results for the ground state energies are obtained. An example 
shows the typical features of the dependence of the ground state energy on the external 
potential. 

Section ^ provides the basis for the rest of this work. However, the sections on the 
determination of heat-kernel coefficients and the calculation of Casimir and ground state 
energies are completely independent and can be read separately. 



2 Developing the formalism: Scalar field on the three dimen- 
sional ball 

To start we focus our interest on the zeta function of the operator (—A + m 2 ) on the three 
dimensional ball B 3 = {x £ IR 3 ; \x\ < a} endowed with Dirichlet boundary conditions. The 
zeta function is formally defined as 

*) = EV. (2-i) 

k 

with the eigenvalues being determined through 

(-A + m 2 )</> fe (x) = \ k 4> k (x) (2.2) 

{k is in general a multiindex here), together with Dirichlet boundary conditions. It is con- 
venient to introduce a spherical coordinate basis, with r = |x| and the angles SI = (0,f)- 
In these coordinates, a complete set of solutions of eq. ( |2.2| ) together with the boundary 
conditions may be given in the form 

(f>l,m, n (r,n) = r~ 1/2 J l+1/2 (wi tn r)Yi m (£l), (2.3) 



7 



with the Bessel function J2+1/2 and Yi m (Q) the spherical surface harmonics JT4|]. The wi :Tl 
(> 0) are determined through the boundary condition by 

Ji+i/2(wi, n a) = 0. (2.4) 

In this notation, using = wf n + m 2 , the zeta function can be given in the form 

00 00 

c(^) = EE(^ + 1 )Kn+^ 2 r, (2.5) 

n=0 1=0 

where w^ n is defined as the n-th root of the 1-th equation (p.4|). Here, the sum over n is 
extended over all possible roots w^ n on the positive real axis, and {21 + 1) is the number of 
independent harmonic polynomials, which defines the degeneracy of each value of I and n in 
three dimensions. 

Eq. ( |2.5[ ) may be written under the form of a contour integral on the complex plane, 

00 r dk 8 

C« = £(» + !) J 3- (k 2 + m 2 )- s ^lnJ l+1/2 (ka), (2.6) 



where the contour 7 runs counterclockwise and must enclose all the solutions of Q2.4| ) on the 
positive real axis, see Fig. [I] (for a similar treatment of the zeta function as a contour integral 
76, 56|). Clearly, (d/dk) In Ji+x/%{ka) = a^/ +1 / 2 (^ a )/^+i/2(^ a )) having simple poles 



sec 



at the solutions of eq. ( p.4| ) and, by the residue theorem, definition ( p.5[ ) is reproduced. This 
representation of the zeta function in terms of a contour integral around some circuit 7 on 
the complex plane, eq. (|2.6|), is the first step of our procedure. 



k-plane 



Figure 1: Contour 7 



As it stands, the representation ( |2.6D is valid for 3f?s > 3/2. However, we are especially 
interested in the properties of £(s) in the range ?Rs < 3/2 and thus, we need to perform the 
analytical continuation to the left. Before considering in detail the /-summation, we will first 
proceed with the /c-integral alone. 
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The first specific idea is to shift the integration contour and place it along the imaginary 
axis. For k — > 0, to leading order, one has the behaviour J u (k) ~ k v /(2 u T{v + 1)) and, in 
order to avoid contributions coming from the origin k = 0, we will consider (with v = 1 + 1/2) 
the expression 



r Ah f) 

C(s) =j — (k 2 + m 2 T 3 ^ In (k-\J u (ka)) , 



(2.7) 



where the additional factor k v in the logarithm does not change the result, for no additional 



pole is enclosed. Using the relations J v {ik) 
one then easily obtains 



e mv J v {-ix) and I v {x) 



-ix-n/2 j 



IX 



C(s) = I dk [k 2 



7T 



m 



Ok 



In (k- u I u (ka)) 



(2.8) 



valid in the strip 1/2 < < 1. Here, the upper restriction, ?R.s < 1, is imposed by the 
behaviour {k — m)~ s of the integrand at the lower integration bound k — » m. For k — > oo one 
uses instead I v {k) ~ e k j\f 7 hxk to find the behaviour k~ 2s and, thus, the restriction 1/2 < ?R,s. 
It gets clear that, with m = 0, this representation is defined for no value of s and a slightly 
more difficult procedure has to be used [|78], ^] . We prefer to include the mass m and consider 
the limit m — > (whenever needed) at a stage of the calculation where the limit will be well 
defined. 

Given that the interesting properties of the zeta function (namely nearly all heat-kernel 
coefficients of (—A), the determinant and the Casimir energy) are encoded at the left of the 
strip 1/2 < < 1, how can we find its analytical continuation to the left? As explained, 
the restriction 1/2 < is a result of the behaviour of the integrand for k — » oo. If we 
subtract this asymptotic behaviour from the integrand in ( |2.8[) , the strip of convergence will 
certainly move to the left, and the hope will be that the asymptotic terms alone can be treated 
analytically, whereas C u ( s ) with the asymptotic terms subtracted can be treated numerically. 
In detail, one needs to make use of the uniform asymptotic expansion of the Bessel function 
I u (k) for v — > oo as z = k/v fixed which, as we will see, guarantees that integration as 
well as summation lead to analytical expressions. One has 



IJuz) 



1 



V2 



TTU (1 + Z + 



k=l 



Uk(t) 



(2- 



with t = 1 / V 1 + z 2 and rj = y/l + z 2 + ln[z/(l + Vl + z 2 )]. The first few coefficients are 



listed in pQ[; higher coefficients are immediate to obtain by using the recursion [80] 



t 

u k+1 (t) = l -t 2 {l - t 2 )u' k (t) + ljdr(l- 5t 2 H(t) 



(2.10) 



starting with uo(t) = 1. As is clear, all the u k {t) are polynomials in t. Furthermore, the 
coefficients D n (t) defined by 



In 



1 ~T~ 2^ „k 



k=l 



Dn{t) 

n=l 



(2.11) 
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are easily found with the help of a simple computer program. For example one has 



D l {t) 



1 5 o 

-t 1 3 , 

8 24 ' 



D 2 (t) = -^--t^-t' 
w 16 8 16 



(2.12) 



By adding and subtracting iV leading terms of the asymptotic expansion, eq. ( 2.11| ), for 

(2.13) 



oo, eq. (2J3) may be split into the following pieces, 



N 



COO = z v {s) + 



with the definitions 

Z v {s) ■ 



sm(7rs 



7T 



dz 



ma/u 



ZV 

a 

-In 



m 



^ {In 



(2.14) 



_V2^(l + z2)|_ 



JV 

£ 

n=l 



£>n(t) 



and 



All 



sin its 



TT 



dz 



ma/ v 



sin(7rs) 



7T 



dz 



ma/ v 



sin(7rs) 



7T 



dz 



ma/ v 



ZV 

a 

ZV 

a 

zv\ 
a J 



m 



m 



m 



d_ 

dz 







In (z- v e vr >) , 



ln(l + z 2 )"4. 



dz 

d_ (DAD 

dz 



(2.15) 
(2.16) 
(2.17) 



The essential idea is conveyed here by the fact that the representation ( 2.13| ) has the 
following anticipated properties. First, by considering the asymptotics of the integrand in 
eq. (2.14) for z — * ma/v and z — > oo, and by considering the behaviour of Z v (s) for v — * oo, 
which is v~ 2s ~ N ~ 1 , it can be seen that the function 



Z( S )=£(2/ + 1)Z^) 

1=0 



(2.18) 



is analytic on the half plane (1 — N)/2 < ?R.s. (The integral alone has poles at s = k S IN 



coming from the z — > ma/v region which is seen by writing [(zv/aY 



m 



2i -s 



(-l)^(T(l 



s)/r(j + 1 — s))(d J /d(m 2 y)[(zv/a) 2 — m 2 ]~ s+: > . But these are cancelled by the zeroes of 
the prefactor.) For this reason, it gives no contribution to the residue of in that range. 
Furthermore, for s = —k, k £ JNq, k < (—1 + N)/2, we have Z{s) = and, thus, it also 
yields no contribution to the values of the zeta function at these points. This result means 
that the heat kernel coefficients, see eqs. ( |1.12 ) and ( 1.13j ), are just determined by the terms 
Ai(s) with 



A(s) = J2(2l + l)AUs). 

1=0 



(2.19) 
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As they stand, the A^(s) in eqs. fl2.15| ), J2. 16 ) and ( p.!7| ) are well defined on the strip 
1/2 < < 1 (at least). We will now show that the analytic continuation in the parameter s 
to the whole of the complex plane, in terms of known functions, can be performed. Keeping 
in mind that Di(t) is a polynomial in t (see the examples in eq. ( |2.12| )), all the A\ (s) are in 



fact hyper geometric functions, which is seen by means of the basic relation [77] 

r(c) 



2 F 1 (a,b;c; z) 



T{b)T(c-b) 



f 1 dt t b -\l-t) c - b - l {l-tz) 
Jo 



(2.20) 



Let us consider first in detail A^^s), Aq(s), and the corresponding A-i(s), Aq(s). One 
finds immediately that 



i n ■-' r ( s 2) 

ram 7^rr~, — — 2F1 



-2..s 



2^F 



-— m 
2 



r 00 



1 

2' 




-(-)' 

2' 2' \maj 



(2.211 



--m- 2s 2 F l I 1, s; 1; 



(-)' 

\maj 



-—m 

4 



-2s 



1+1- 

ma 



(2.22) 



where in the last equality we have used that 2F1 (b, s; b; x) = (1 — x)~ s . These representations 
show the meromorphic structure of A l ^_ l (s), Aq{s) for all values of s. 

The next step is to consider the summation over /. For A^^s) this is best done using a 
Mellin-Barnes type integral representation of the hypergeometric functions, namely 



2 F l (a, b; c;z) 



r(c) 



1 



T(a)T(b) 2m 



dt 



T(a + t)T(b + t)T(-t) t 
T(c + t) 1 Zj ' 



(2.23) 



where the contour is such that the poles of T(a + t)T(b + t)/T(c + t) lie to the left of it and 
the poles of r(— t) to the right [77]. The contour involved in eq. (|2.21| ) is shown in the Fig. ^. 

It is clear that one wishes to interchange the summation over / and the integration in 
(|2.23j) in order to arrive at a Hurwitz zeta function, which is defined as 



1=0 



> 1. 



(2.24) 



However, as is well known, one has to be very careful with this kind of manipulations, 

]. Before interchanging 



what has been realized and explained with great detail in [81, 32 



the summation and integration one has to ensure that the resulting sum will be absolutely 
convergent along the contour C. Applying this criterium to A_i(s), 



A-i(s) 



£(21 + 1) 

1=0 



m 



-2s 



2v^ 



am 



r(.-i) 
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1F1 



1 
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'l + l 



ma 



2 -m- 2s 
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t-plane 



Figure 2: Contour C for eq. (2.21) 



it turns out that we may interchange the J2i an d the integral in eq. ( g2jp only if for the 
real part of the contour the condition 5RC < — 1 is satisfied. However, the integrand 
r(-l/2 + t)T(s - 1/2 + i)/r(l/2 + i) has a pole at t = 1/2 (assume for the moment 5Rs > 1 
so that all arguments go through). Thus the contour C coming from — ioo must cross the 
real axis to the right of t = 1/2, and then once more between and 1/2 (in order that 
the pole t = of r(— t) lies to the right of it), before going to +ioo (see Fig. 2). That is, 
before interchanging the sum and the integral we have to shift the contour C over the pole 
at t = 1/2 to the left, which cancels the (potentially divergent) second piece in A_i(s). This 
term, —(I + l/2)m~ 2s /2, is the result of the factor lnk~ u introduced in eq. (2.7) to avoid 



contributions coming from the origin, and its crucial importance is made explicit in the above 
step. 

Closing then the contour to the left, we end up with the following expression in terms 
of Hurwitz zeta functions (the contour at infinity does not contribute, which is seen by 
considering the asymptotics of the integrand) 

2s oo p (j + s - ±) 

For Aq(s), one only needs to use the binomial expansion in order to find 
2s «3 f_iV 

MS) = ~2T(s) £ { —f(rna)VT(s + j)( H (2j + 2,-1; 1/2). (2.26) 

The series are convergent for \ma\ < 1/2. Finally, we need to obtain analytic expressions 
for the Ai(s), i G IN. As is easy to see, they are similar to the ones for A-i(s) and Aq(s) 
above. We need to recall only that Di(t), eq. (2.11), is a polynomial in t, 



D t {t)=Y J ^bt l+2 \ (2.27) 

6=0 
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which coefficients are easily found by using eqs. (2.10) and (2.11) directly. Thus the 
calculation of A"(s) is essentially solved through the identity 



dz 



ma/ v 



ZV 

a 



m 



— 

dz 



-m 



-2s n r( s + §)r(i- s ) 

2(mo)" r (1 + § ) 



1+1- 

ma 



(2.28) 



The remaining sum may be done as mentioned for Aq(s), and we end up with 



Ms) 



2a l 



T(s) 



00 f— IV 



( H (-l + i + 2j + 2s;l/2) 



3=0 



6=0 



r s + 6 + j + 



r[6+§ 



(2.29) 



convergent once more for \ma\ < 1/2. Restricting attention to the massless field, the asymp- 
totic contributions take the surprisingly simple form 



,2s 



A (s) 
Ms) 



F(s - 1/2) 

20F r(s + l) 



Ch(2s-2;1/2), 



,2s 



2 

2a 2s 



£h(2s - 1; 1/2), 



(2.30) 



— -Cif(i-l + 2 S ;l/2)^^, 6 



r (s + 6 + j + 
r(^ + §) 



In summary, the analytical structure of the zeta function for the problem considered is 
made completely explicit. Eq. ( [2.30 ) allows for the direct determination of residues, function 
values and derivatives at whatever values of s needed. Clearly, all ingredients can be easily 
dealt with by Mathematica and calculations can be easily automized. The remaining piece 
to obtain the full zeta function, eq. ( p. 18 ) and ( |2.14| ), is suitable for numerical evaluation 
for values of s in the half plane mentioned. Furthermore, an analytical treatment of Z'(0) 
is possible, which leads to the determinant C'(0) °f t ne Laplacian on the ball in terms of 
elementary Hurwitz zeta functions 84, |5j|. For various different approaches in a similar 
context see also 36, [7(1, S^, j§8|, with possible applications in quantum cosmology |8^, 9C]. 



3 Calculation of heat-kernel coefficients via special cases 

Let us now come to the first application of our analysis in Section B namely the calcula- 
tion of heat-kernel coefficients for Laplace-like operators on smooth manifolds with smooth 
boundaries. Our main emphasis is on the determination of the boundary contribution to the 
heat kernel coefficients. This is motivated because it is here that the special case of the ball 
gives rich information. But it is also justified because the calculation of the volume part 
is nowadays nearly automatic [91, p2l |93| and these terms do not depend on the boundary 



conditions [94] and are thus known already for all problems to be dealt with. 
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Concerning the relevance of heat-kernel coefficients, in the Introduction we have already 
mentioned that a knowledge of these coefficients is equivalent to a knowledge of the one- 
loop renormalization equations in various theories p^] , which provides one reason for the 
consideration of heat-kernel coefficients in physics. In addition, if an exact evaluation of 
relevant quantities is not possible, asymptotic expansions (with respect to inverse masses, 
slowly varying background fields, high temperature,...) are often very useful and naturally 
given in terms of heat-kernel coefficients. An especially important link between physics and 
mathematics is provided by index theorems ||, again, with the well-known connection to 
the heat equation. But, in mathematics, the interest extends to basically all of Geometric 
Analysis, including analytic torsion [|l^, 11], characteristic classes j9j, sharp inequalities of 



borderline Sobolev and Moser-Trudinger type 95], etc. 

We will start the Section by giving what can be named the general form of the heat-kernel 
coefficients for Dirichlet and Robin boundary conditions. (We will explain the approach for 
these boundary conditions; comments on various other boundary conditions are given at the 
end of the Section.) As we will see these are built from certain geometrical invariants with 
unknown numerical coefficients. Relations between the unknown coefficients can be derived 
by conformal transformation techniques most systematically used by Branson and Gilkey 96 ] . 
However, in order to determine the numerical coefficients, additional information is needed. 
The product formula gives a certain subset of the numerical coefficients but in general by far 
not enough to complete the calculation by using the conformal techniques |9£|| . Especially 
the group of terms containing the extrinsic curvature is not even touched by the product 
formula and our calculation on the ball will turn out to be very valuable. Combined with 
the conformal techniques, the application of index theorems, and additional examples, the 
conglomerate of all methods allows for the determination of (at least) the leading heat-kernel 
coefficients for all classical boundary conditions. 

Let us clearly state that the coefficients we are going to "determine", namely the coeffi- 
cients up to 03/2 for Dirichlet and Robin boundary conditions, are known already for a long 
time, see. f.e. |9(|. It is just to explain the ideas in detail that we have chosen the first coef- 
ficients and the simplest boundary conditions. But in fact, the method has been successfully 
applied to containing more than 100 terms |97j (see also |98|| ), showing the effectiveness 
of the procedure employed. 



3.1 General form of the coefficients for Dirichlet and Robin boundary con- 
ditions 

In order to be as self-contained as possible, in the following we are going to summarize some 
basic properties of heat-kernel coefficients on manifolds with boundary. We follow p6| ]. 

Let A4 be a compact D = {d + l)-dimensional Riemannian manifold with boundary dAi. 
Let V be a smooth vector bundle over M. equipped with a connection and finally let E 
be an endomorphism of V. Our interest is then in Laplace- type operators of the form ( |1.4[) , 

P = - 9 ij vYvJ - E. (3.1) 

Let us mention that every second order elliptic differential operator on M with leading 
symbol given by the metric can be put in this form. We will see this explicitly below, starting 
with eq. (|I|). 

If there is no boundary, it is well known that P defines a symmetric operator. If, however, 
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there is a boundary present, Green's theorem says 

(v, Pw) L 2 — (Pv, w) L 2 = I dx(v^ Pw — (Pv^w) = / dy(v J m w — v^w- m ), (3.2) 

JM JdM 

with dx and dy the volume elements on A4 and dM, and v ]m the normal covariant derivative 

of v with respect to the exterior normal N to the boundary dj\A. In order that P be symmetric 

one has to impose boundary conditions. Obvious possibilities are the classical Dirichlet or 

Robin boundary conditions, 

B~4> = (j)\ dM and B^cj) = (4> ;m - S(p) \ dM , (3.3) 

with 5 a hermitian endomorphism of V defined on dj\A . Clearly, in these cases the argument 
of the integral over the boundary in (|3,2j ) vanishes and, with these boundary conditions, a 
symmetric operator P is defined. 

Another possibility is to impose a mixture of Dirichlet and Robin boundary conditions. 
In order to do this, a suitable splitting of V = V- €B is needed and in V- one imposes 
Dirichlet and in V+ Robin boundary conditions. Also here, the integrand itself vanishes. 

However, a further possibility is that the integrand in (3.2) does not vanish but equals a 



boundary divergence. This condition involves tangential derivatives and, in the mathematical 
literature, it is sometimes referred to as oblique. Although these last boundary conditions 
have been the subject of classical analysis (see, f.e. |99|, IOC, |101| ), very little is known 



about the associated heat equation asymptotics. We will comment on the application of the 
techniques to all these boundary conditions at the end of the section, but concentrate now 
on conditions (^). 

To have a uniform notation, we set 5 = for Dirichlet boundary conditions and write 
Bf. If F is a smooth function on A4, there is an asymptotic series as t — > of the form 



Tr L2 (j 



Fe- tp ) «X>-fa n (f,i>B£), (3.4) 

n>0 



where the a n (F, P,Bg) are locally computable Q. 

The use of the smearing function F is important for at least three reasons. First, near 
the boundary the heat-kernel behaves like a distribution and, by studying Tr i 2 (Fe~ tp ), this 
local behaviour is recovered. As an example consider <9A4 = 0. For F = 1 volume divergences 
are integrated away, which is not possible if a smearing function is present. A second reason 



is that for the functorial formalism to be described [96], the smearing function is at the heart 



of this method. Finally, it is exactly this smeared coefficient appearing in the integration of 



conformal anomalies that is relevant for several physical applications [102, 103, 104, 105 , 106 , 



i07| , pi |i09 , |iiq , |m| l 



To state the general form of the coefficients in ( |3.4[ ), let us introduce some notation. Here 
and in the following F[M] = Tr J M dx F(x) and F[dM] = Tr J dM dyF{y), with Tr the fiber 
trace. In addition, ";" denotes differentiation with respect to the Levi-Civita connection of 
j\A and ":" covariant differentiation tangentially with respect to the Levi-Civita connection 
of the boundary. Furthermore, Q is the curvature of the connection V^, [Vj^Vj] = O^, 
and Rijki, R%j, R, are as usual Riemann tensor, Ricci tensor and Riemann scalar. Finally 
let N U (F) = F- m ,„ be the u th normal covariant derivative. Then there exist local formulae 



a n (x,P) and a n ^{y,P,B^) so that [112] 



2n-l 



a n (F,P,BT) = {Fa n (x,P)}[M} + { £ N"(F)a n>v (y, P, B*)}[dM]. (3.5) 

u=0 
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Important homogeneity properties follow from the Seeley calculus so one has for 
< c e IR Iffli 



a n (x, c~ 2 P) = c- 2n a n (x, P), a n , v {y, c~ 2 P, c~ x B^) = c^ 2 ^ a n>v {y , P, Bj). (3.6) 

Physicists would say that P carries dimension length -2 , thus in order that e~ tp makes 
sense, t has to have dimension of length 2 . Then, eq. is dimensionless and a n (F,P,B^) 
must have dimension of length~ D+2n . As a result, a n (x,Pis) has dimension of length to the 
power 2n and a Ujl/ (y, Pb-i&s) to the power 2n — v which is equivalent to the above. 

The interior invariants a n (x,P) are built universally and polynomially from the metric 
tensor, its inverse, and the covariant derivatives of R,£l, and E. By Weyl's work on the 
invariants of the orthogonal group 1 1 13 ] , these polynomials can be formed using only tensor 



products and contraction of tensor arguments (indices). If A is a monomial term of a n (x, P) 
of degree (&r, kn, fee) in (R, ft, E), and if ky explicit covariant derivatives appear in A, then 
by the homogeneity property of a n (x, P), 

2(k R + k n + k E ) + k v = In. 

When considering the boundary invariants a n ^{y,P,B^) we must also introduce the sec- 
ond fundamental form K ab = (V ea eb, N), K = K%, where {ei,...,e<i} is an orthonormal 
frame of T(dM), the tangent bundle of the boundary and, when considering Robin bound- 
ary conditions, we must also consider the tensor S. Given that these are defined only at 



the boundary, we only differentiate {K, S} tangentially. We use Weyl's [113] theorem again 
to construct invariants. The structure group now is 0(D — 1), and the normal N plays a 
distinguished role. If A is a monomial term of a njV (y, P,B^) of degree (kn,kn,kE,kK,ks) 
in (i?, O, E, K, S), and if k\? explicit covariant derivatives appear in A, then once more by 
homogeneity 

2(k R + kn + k E ) + k K + k s + ky = 2n - v. 

By constructing a basis for the space of invariants of a given homogeneity, we write down 
the following general form of the heat-kernel coefficients 



a (F,P,B%) = {^r D ' 2 F[M], (3.7) 
a 1/2 (F,P,B%) = 5(47r)- d / 2 F[dM], (3.8) 
ai(F,P,B$) = (Att)- d / 2 6- 1 {(6FE + FR)[M] + (b FK + b 1 F m + b 2 FS)[dM}}, (3.9) 

a 3/2 (F, P, B|) = ——^ {F{coE + Cl R+ c 2 R mm + c 3 K 2 + c 4 K ab K ab + c 7 SK+ csS 2 ) 

+ F. m (c 5 K + c 9 S) + ceF mm }[dM]. (3.10) 

All numerical constants involved have the very important property of being independent 



of the dimension D [112]. A direct proof has been given in 96] and this goes as follows. 
First one shows a product formula for heat-kernel coefficients. Let Ai = «Mi x Ai 2 and 
P = Pi®l + l(£iP 2 and dM 2 = and let S only depend on coordinates in A4\. Then 

a n (x,P) = a p (x 1 ,P 1 )a q (x 2 ,P 2 ), 

p+q=n 

anAv^ p ^s) = a Pi „(y 1 ,P 1 ,Bj)a q (x 2 ,P 2 ). (3.11) 

p+q=n 
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This is a purely formal computation because by separation of variables the heat-kernel of 
the operator P gets the product of the heat-kernels of Pi and P2 and the result follows by 
just comparing powers of t. 

To avoid the appearance of factors of \^4tt normalize for the moment ao(x,P) = 1. 
Application of formula to (M 2 ,P 2 ) = {S 1 ,-d 2 /d9 2 ) leads to a n (xi,P\) = a n ((xi,9), 



P) and a niU (yi,Pi,B^) = a n<u ((yi, 9), P, Bj), because 00(0,^2) = 1 and a q (9,P 2 ) = for 
q > 0. However, invariants formed by contractions of indices are restricted from Aii x S 1 to 
Ail by restricting the range of summation, but have the same appearance. This shows that 
the numerical constants are independent of the dimension. 

We are thus left with the task of the determination of the universal numerical constants 
by whatever method. An obvious (rich) source of information are special case calculations. 
Let us discuss this in detail for the example of the ball in order to motivate the calculation of 
the coefficients for this setting. The way to deal with the three dimensional ball and Dirichlet 



boundary conditions has been shown in Section g. In Section 3.2 we will show how these 
results can be generalized to arbitrary dimensions D = d + 1. In addition, the way to deal 
with Robin boundary conditions with a constant endomorphism S is explained. For these 
manifolds we will have K h a = 5\. As a result we get K = d, K^K^ = d, K 2 = d 2 , and so 
on. The polynomials, traces and contractions of K a ^ give a polynomial in the dimension d. 
In addition, for the example of the ball, we have R a tce = 0) and we have chosen P = — A_/vi, 
thus E = 0. Finally, we included no smearing function and have F = 1. In this setting, 
restricting Q3.8| ) to the ball, we have 

a 1/2 (l,-A M ,BT) = 5(47r)- d / 2 \S d \, 

with the volume \S d \ = 2ir( d+1 ^ 2 /T((d + l)/2) of the d-spheie. Calculating the heat-kernel 
coefficients explicitly one will find the "unknown" numerical constants 5 and then knows 



ai/ 2 (F, P,Bs) for an arbitrary manifold by eq. (3J3). Passing on to a±, on the ball we have 



ai (l,-A M ,B$) = (4ny D / 2 6- 1 \S d \(b d + b 2 S). 

Just by comparing powers of d and S one can determine 60 and b 2 from the explicit 
ai(l, — A_/k Bg) on the ball. Let us stress that, if we include a smearing function F(r) 
into the formalism, then F ;m = (d/dr)F(r) is the derivative with respect to the exterior 
derivative, and we would have determined also 61. By application of ( |3.9| ) we see that one 
can get a±(F, P,B^) by just having the result on the ball. 

Continuing with the same argumentation, in eq. ( |3.1C| ) one can determine 03,04, C7 and eg 
for F = 1 and furthermore C5, eg and cq including an F(r). Let us stress that C3 and C4 both 
can be determined only because we will perform our calculation in arbitrary dimension (this 
observation gets more important for the higher coefficients). Thus only 3 of 10 unknowns 
are left and it gets clear that the special case calculation chosen contains rich information. 
However, one also realizes, and this was, of course, clear from the beginning that the example 
can not determine the full coefficients, because E and the Riemann tensor vanishes. Both 
aspects can be improved a bit by including a mass (as we have done) and by dealing with a 
so-called bounded generalized cone J59[ |. But the information obtained thereby is also very 
easily obtained by an application of the product formula ( |3.11 ) and we will not give details 



of this generalization here. However, product manifolds have vanishing normal components 
of the Riemann tensor (their appearance starts with a 3 / 2 ) and the corresponding universal 
constants have to be determined by different means. 
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An extremely effective method to continue at this point is to study the functorial proper- 

]. We summarize here the main points in order to 



ties of the heat-kernel coefficients [ 114 



allow for a self-consistent reading. Consider the one-parameter family of differential operators 

P(e) = e~ 2eF P (3.12) 



and boundary operators 



(3.13) 



Here, e is a real- valued parameter and F, as before, is a function. Eq. ( 3.13[ ) guarantees 
that the boundary condition remains invariant along the one-parameter family of operators 
(see below, eq. ( p. 19 )). Then one might ask what is the dependence of the heat-kernel coeffi- 
cients on the parameter e. The relevant transformation behaviour is described by the following 



Lemma: 



(") j e \ e=0 a n (l,P(e),B^(e)) = (D-2n)a n (F,P,B^), 



(b) If D = 2n + 2, then ^ | e=0 a n (e" 2e/ F, P(e), 0|(e)) = 0. 



(3.14) 
(3.15) 



Proceeding formally, (a) is proven by considering 



d 



— L-o Tr L2 



-£Tr 



-tP 



2tTr L 2 [FPe 



-tP 



(Fe-' p 



and comparing powers of t in the asymptotic expansion. For the necessary justification of the 
analytic steps see [112]. Part (b) is much the same by starting with (d/de) | e= oTr i 2(e _2e ^ Fe~ tp ^ e 



This Lemma will be applied in Section |3.3| ; let us here only explain the way it determines 
universal constants. First, we obviously need the heat-kernel coefficients for the operator 
-P(e), so let us see how the family ( 3.12[ ) of operators can be generated. Let P be an arbitrary 
second order differential operator with leading symbol given by the metric tensor. In local 
coordinates we have 



d 2 



dx' l dxi 



+ P k 



d 
dx k 



Obviously, P(e) is obtained by defining 

g ij (e) = e- 2eF g i \ P k (e) 



e -2eFpk 



+ Q 



(3.16) 



-2eF 



Q. 



In order to obtain the heat-kernel coefficients of P{e) in the form (3.7) — ( 3. lCj ) , write P 



invariantly in the form (|3.1[), P = — (<7 y Vjv*f + E). Let lo\ be the connection 1 form and 



Comparing the two different representations of P one finds 



consequently fli 



U>i t j + U>iUJj 



ujjUJi with "," 



E 



l - 9ll (P* + g^T% ) . 



Q-g lJ ( 



'■j 



J 



the partial derivative. 

(3.17) 
(3.18) 
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with the Christoffel symbols r*- fc . As a result, this defines 

w,(e) = ui + ^e(2 - D)F, l (3.19) 

fiy(e) = fiy (3.20) 

£(e) = e- 2eF (^ + i( J D-2) e A A4 F+i(D-2) 2 e 2 F. fe F. fc ). (3.21) 



This shows that the leading a n (f,P(e),B^(e)) are given by eqs. (|3.7l) — ( 0.10 ) once the 



above definitions are used and once all geometrical tensors and covariant derivatives are 
calculated with respect to the metric gij(e) (for many useful relations see |)(J). Let us here 
only mention that Dirichlet boundary conditions are obviously conformally invariant, and 
that with 

S(e) = e-* ( S - e^F, 



the same holds for Robin conditions. This is seen immediately from eq. ( 0.19 ) and the bound- 



ary condition (|Q|). As a result, (d/de)\ e =oak(l,P(e),Bg(e)) will have the same appearance 
as ak(F,P,Bg), and Lemma ( 3.14[) will give relations among the universal constants, as well 



as ( |3.15|) does. To make it very clear, look a t a 3 / 2 (F, P, Bg). Then (d/de)E(e) contains a 
term F- mm being part of AF, see ( 3.21 ). Eq. ( 3.14| ) then states that the numerical constant 



Co is connected with cq (more invariants are involved however). Let us stress already here, 
that due care must be taken that only independent terms are compared in eqs. ( p. 14 ) and 



( 3.15| ) and that partial integrations (or more involved manipuliations) may be necessary to 



see that apparently independent terms are actually dependent. 

On its own, the functorial method is unable to determine the coefficients fully. But, given 
a subset of numerical coefficients, found f.e. by special case calculations, the method provides 
the required information with relative ease. For that reason, we start by applying the product 
Lemma ( 3. 11[ ) and, by calculating the heat-kernel coefficients for the Laplacian on the ball, 



determine (a subset of the) universal constants and complete the calculation by use of the 
functorial properties. 

In the last step it will turn out that the generalization of the calculations in Section 
3.2 to the smeared zeta function is essential. This is, as seen in eq. ( 3.14| ), because the 



functorial techniques (apart from other things) yield relations between the smeared and non- 
smeared case; the information one can get on the "smeared side" is crucial to find the full 
"non-smeared" side. The way this can actually be done is explained in [pi]] . 

3.2 Scalar field on the .D-dimensional ball 

As briefly explained, for the applications of our calculations to the heat-equation asymptotics, 
it will be very important that results in arbitrary dimensions are available. For that reason, 
let us briefly explain how the three dimensional calculation of Section ^ is generalized to 
arbitrary dimension D = d + 1 J5J|. We put the radius of the ball a = 1, the dependence of 
the results on the radius are easily recovered by dimensional arguments. 

Instead of (^^) , the nonzero eigenmodes of A that are finite at the origin have eigenvalues 
—a 2 and are of the form 

^j2 Y l + D/2(n), (3.22) 
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d-l 



2 



where the spherical harmonics satisfy 

&nYi + d/2(V) = -XfY l+D/2 (n) (3.23) 
with Ajv - the Laplacian on the sphere, and 

v 2 = \f + (d- l) 2 /4. (3.24) 

The eigenvalues are [74] 

Af = /(/ + d-i) = (/ + ^) 2 -(^-) 2 , leJNo, 

such that 

v 2 = (/ + 
and they have degeneracy 

Let us next impose the boundary conditions and this time let us consider Dirichlet as well 
as generalized Neumann (or Robin) boundary conditions. In the notation of, for example, 
[115], these read explicitly 

J u (a) = (3.25) 

for Dirichlet and 

uJ v (a) + aJ' v (a) = (3.26) 

for Robin, with u 6 IR. In the following we will write u = 1 — y — (3, such that (3 = cor- 
responds to Neumann boundary conditions and (3 is to be identified with the endomorphism 
S. 

As we have seen in the three dimensional calculation, the angular momentum sum leads 
to a zeta function of the type, 

Cat(s)=£^K 2S > (3.27) 

and we anticipate this to be the central object to state the asymptotic contributions. Its 
definition is clearly motivated by the calculation in 3 dimensions where 2£#(2s — 1; 1/2) 
plays the role of Cm( s )- 

Our first aim will be to express the whole zeta function on the ball 



cw = E 



a~ 2s , 



as far as possible in terms of this quantity. That is, we seek to replace analysis on the ball 
by that on the sphere. 
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Following the analysis of the previous section for Dirichlet boundary conditions, the start- 
ing point is the representation of the zeta function in terms of a contour integral 



(3.28) 



where the anticlockwise contour 7 must enclose all the solutions of ( |3.25 ) on the positive real 



axis. 



As we have already seen, it is very useful to split the zeta function into two parts, 



N 



C(s) = Z(s)+ £ m*)- 



(3.29) 



i=-l 



Performing identical steps as before, the different pieces are determined to be 

00 

z{s) = *j^Y,d(v) j dz(z V )- 2s ^ z (\n(z-»I y {zv)) 



(3.30) 



In 



,,!/// 



and 



A (s) 

M*) 



r s 



40F r( s + i) 



.V2^(i + Z 2)i 



Cat (s- 1/2), 



N 

E 

71=1 



D n {t) 



\qm{s) 



1 \ s > b=0 1 



r ( s + b + 1/2) 

(6 + V2) ' 



(3.31) 
(3.32) 
(3.33) 



These results can be read off from eqs. ( |2.14| ) and ( |2.18|) once the definition ( |3.24j ) for u 2 
is used and 2£#(2s — 1; 1/2) is replaced by Ov(s), eq. (|3.27|) . The function Z(s) is analytic 
on the strip (d — 1 — N)/2 < fits, which may be seen by considering the asymptotics of the 
integrand in eq. ( 3.30| ). 

As is clearly apparent in eq. (^3l|)-(p3D, sph ere contributions are separated from radial 
ones. Again, the analytical structure is made very explicit and the result is very well organized 
by the introduction of the zeta function 0v( s )- 

In order to treat Robin boundary conditions, only a few changes are necessary. In addition 
to expansion (2J3) we need [|116| , [SCf] 
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i + E 
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(3.34) 



with the Vk{t) determined by 

v k (t) =u k (t)+t{t 2 -I) 



:Uk-i(t) + *«fc_i(i) 



(3.35) 
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The relevant polynomials analogous to the D n (t), eq. (|2.11|), are defined by 



In 



v k (t) , 1-D/2-P 



+ 



Ufc(i) 



fc=l 



fc=l 
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n=l 



Af n (t) 



(3.36) 



and have the same structure, 



M n (t) =Y,Zn,bt 



n+2b 



(3.37) 



One may again introduce a split as in eq. ( |3.29| ) with A^(s) = A^i(s) and Aq(s) = 
— Aq(s), where the upper index R indicates that these are the results for Robin boundary 
conditions. The A^(s) are given by eq. ( 3.33| ) once the Xi ja is replaced by zi^ a . In addition 
one finds 

oo 

z R {s) = s ^^Y,d(u) J dz{zv)- 2s ^\n{{l-D/2-(3)I u {zv)+zvIl(zv)) 



In 
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M n {t) 



(3.38) 



The remaining task is the analysis of the meromorphic structure of the zeta function 
Cv(s). It reads explicitly 



CtfW = £ (a+- _ 1) C£-ra!(i + ''- lx " 2 " 
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l\(d-l)\ 



(3.39) 



Writing 



d(l) 



l + d-1 
d-1 



+ 



l + d-2 
d-l 



it is seen immediately that this is a sum of Barnes zeta functions [117, 118, 119] defined as 

1 



Cb(s,o) 



m=0 



(a + mi + ... + m d ) 
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(3.40) 



In detail, one finds 



(Ms) = Cb 2s, 



d + 1 



+ Cb 2s, 



d-l 



(3.41) 



Its residues are determined in the following way. Using the integral representation 
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2^ 



where L is the Hankel contour, one immediately finds for the base function 
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For the residues this yields (m = 1, 2, d) 



d D {d-i) 
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d—m 



(d - \){m - 2)\{d - m)V 



with the D {d 1] defined through (cf pOR ) 
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(3.44) 



(3.45) 



Obviously D 



(d-l) 



for v odd, so there are actually poles only for m = 1,2, with 



d — m even. The advantage of this approach is that known recursion formulas allow for an 
efficient evaluation of the as polynomials in d, [121]. 

Using eq. ([3.44 ) in eqs. (|3.31 ) — 03.33|) we find, by the connection ( 1.12j ) for the heat-kernel 
coefficients a k j 2 with Dirichlet boundary conditions, 
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(3.46) 



where (y) n = T(y + n)/T(y) is the Pochhammer symbol. Eq,( |3.46 ) exhibits the heat-kernel 
coefficients as explicit functions of the dimension d (partly encoded in ) . Clearly, evalu- 
ation of eq. ( |3.46| ) is a simple routine machine matter, because all ingredients can be found 
by simple algebraic computer programs. 

For later use, the polynomials up to 0,3/2 are listed in the following, 

(4vr) d / 2 1 
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(4vr) D / 2 



1/2 



384 



(3.47) 
(3.48) 
(3.49) 



\S d \ 
(4vr) d / 2 

-|^|- a 3/2 

For Robin boundary conditions one has to make the modifications outlined above eq. ( |3.34| ) . 
The results, up to a 3 / 2 are listed below, 
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(3.50) 
(3.51) 
(3.52) 
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Further coefficients could be calculated with ease, f.e. for the first 20 coefficients the program 
needs about 2 minutes. 



3.3 Determination of the general heat-kernel coefficients 

We now compare, one by one, the general form of the coefficients with our special case 
evaluation. The coefficient ao is, by normalization, 

a (F,P,BT) = (4n)- D / 2 F[M}. 

The next one is 

a 1/2 (F,P,B$) = S(A7r)- d/2 F[dM}. 

For the ball this means 

a 1/2 (F,-A M ,B^) = 5(47r)- d / 2 F(l)\S d \. 
Using the relations ( |3.47| ) and ( |3.50| ) we can immediately determine 8, 

The coefficient a\i 2 is thus given for a general manifold from the result on the ball (which 
was clear of course). Passing on to a\, the general form is 

oi(F,P,Bj) = (ATT)- D / 2 6- 1 {(QFE + FR)[M] + (b FK + b 1 F. m + b 2 FS)[dM]} 

In our special case on the ball, K\ = 8\ and thus 

ai (F, -A M ,g$) = (4ir)- D / 2 6- 1 vol(S d ) {b F(l)d + hF, r (l) + b 2 F(l)S} . 

Comparing with the results ( |3.48| ), ( |3.51 ), one finds 

b = 2, b 2 = 12. 



Taking into account also the smeared calculation [61], in addition one can find 

h = (3- -3+). 

Thus, our special case also gives the entire a\ coefficient without any further information 
being needed. It is very important that the calculation can be performed for an arbitrary 
ball dimension D, and also for a smearing function F(r). This allows one to just compare 
polynomials in d with the associated extrinsic curvature terms in the general expression and 
simply to read off the universal constants in this expression. 

We continue with the next higher coefficient, with the general form, 



a 3/2 (F,P,Bf) = 96(47r)d/2 [F(c E + Cl R+ c 2 R mm + c 3 K 2 + c 4 K ab K ab c 7 SK+ csS 2 ) 

+ F. m (c 5 K + c 9 S) + c 6 F. 

in 
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The (smeared) ball calculation immediately gives 7 of the 10 unknowns, 



C3 = (7-,13 + ), c 4 = (-10/,2+), c5 = (30-,-6 + ), 
c 6 = 24, c 7 = 96, c 8 = 192, c 9 = -96. 



Next, apply the lemma on product manifolds ( 3.11 ) [pq ]. For a 3/ / 2 this means 
as/2{y,P,B^) = a 3 / 2 (yi, Pi,B^)a (x 2 , P 2 ) + a 1/2 (y!, P 1 ,B^)a 1 (x 2 , P 2 ). 
We will choose P\ = — Ai and P2 = — A2 + E(x 2 ), with obvious notation, to obtain 

Sm^icoE + c 1 R(M 2 )) = 66-\eE + R{M 2 )), 
where we used in addition R(M\ x M2) = R(M\) + R(M 2 ). This gives 

Co = 96, ci = 16. 

It is seen that the determination of 03/2 is relatively simple, once the ball result is at 
hand. The lemma on product manifolds is also very easily applied and already only one of 
the universal constants Cj, namely c 2 , is missing. 

The remaining information is obtained using the relations between the heat- kernel coeffi- 
cients under conformal rescaling, fl3.14|) . Setting to zero the coefficients of all terms in ( [3.14| ) 
gives several relations between the universal constants Cj. We will need only one of them. 
Thus, setting to zero the coefficient of F- mm gives 

I(D _ 2) Co - 2(D - l)a — (D — 1)02 ~(D- 3)c 6 = 

and so C2 = —8 for Dirichlet and Robin boundary conditions. This completes the calculation 
of a 3/2 . 

It gets clear, already, that the combination of the different methods is extremely effective 
in obtaining heat-kernel coefficients and, as mentioned, a 5 / 2 has been obtained proceeding 
in this way. For other boundary conditions, the general form of the heat-kernel coefficients 
can contain several other geometrical quantities, as projectors for mixed boundary conditions 
(for a discussion of their physical relevance see |122(| ) or tangential vector fields for oblique 
boundary conditions. Although the analysis gets partly more involved, the ideas described 
seem to apply generally and, as mentioned, many new results have already been obtained 
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4 Casimir and ground state energies 



In this section we apply the techniques developed in Section g to the calculation of vacuum 
energies when spherically symmetric boundaries or external potentials are present. The 
boundary calculation is an immediate continuation of Section [2]. When external potentials 



are present, the implicit eigenvalue equation (2.4) is replaced by an asymptotic equivalent, 
see eq. ( 4.32 ), which allows the energy to be expressed by quantum mechanical scattering 
dates. 
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4.1 Massive scalar field with Dirichlet boundary conditions 

So let us start to analyse the influence the mass of a field has on the Casimir energy. We 
have seen that the Casimir energy is generally plagued by ambiguities. For a massive field, 
however, there is the possibility to separate (at least in principle) the classical part and the 
quantum part of the energy. The idea is that, if the mass of the quantum field tends to 
infinity, quantum fluctuations should die out and as a result, in this limit, the quantum 
contributions to the Casimir energy should vanish. As we will see, this provides a unique 
definition 126]. Before we proceed, let us first discuss this in more detail. 

To actually impose the condition, the m — > oo behaviour of E vac is needed. This is easily 
derived using the heat-kernel expansion. For the scalar field, the heat-kernel of (—A + m?) 
clearly is 



K(t) 



with K-A,(t) the heat-kernel of minus the Laplacian on the ball. The asymptotic expansion 



for m — > oo is then found by employing eq. (1.9) 
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(4.1) 



About a = —1/2 it reads, 
C(-l/2 + s) = 

m 2 ( 1 
+—=a 1 1 + In 

2y/TT \ s 

1 / 1 

~7T~f =a 2 --2 + ln 

Thus, in order to impose the normalization condition 




+ C(l/m) +(D{s). 



m — »oo 



o. 



(4.2) 
(4.3) 

(4.4) 



this defines the terms to be subtracted. As is seen, the zeta functional regularization used 
leaves the contributions of the coefficients with half integer index finite in the limit s — > 0. 
This is a specific feature of this regularization. However, in other regularizations, as for 



example the proper time cutoff [28| or the exponential cutoff 1 127] , these contributions are 
divergent when the cutoff is removed. For this reason, and in order to impose the normaliza- 
tion condition (|4.4j ), we include them among the "divergent" terms suffering renormalization. 
By this, we are led to the following definition of the renormalized Casimir energy, 



with 
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To interpret the subtraction in (4.5) as a renormalization of "bare" parameters, we need 
to consider the following physical system composed of two parts: 

1. A classical system consisting of a spherical surface of radius a. Its energy reads: 

h 

Edass =pV + aS + Fa + k + -, (4.7) 

a 

where V = f Tra 3 and S = 4-Ka 2 are the volume and surface, respectively. This energy 
is determined by the parameters p = pressure, a = surface tension, and F, k, and h, 
which do not have special names. 

2. A quantized field <p(x) whose classical counterpart obeys the Klein-Gordon equation 

(□ + m 2 )<p(x) = 0, (4.8) 

as well as suitable boundary conditions on the surface ensuring selfadjointness of the 
corresponding elliptic operator on perturbations. We choose Dirichlet boundary condi- 
tions as the easiest to handle. 

For this system one can consider three models, which will behave in a different way [126]. 
These models consist of the classical part given by the surface and 

(i) the quantized field in the interior of the surface, 

(ii) the quantized field in the exterior of the surface, 

(iii) the quantized field in both regions together, 

respectively. We will give here the analysis for model (i); models (ii) and (iii) can be treated 
along the same lines [126] and only a few comments will be added. 



The heat-kernel coefficients needed to impose the normalization (4.4) are well known (see 
for instance [12S], or, alternatively, the calculation provided in Section 3.2), and read 

1 3 1 2 1 

ao = ^«, a 1/2 = --a, a x = ^=a, 

"3/a = -24> a2 = ^7^- 

Thus, we have five divergent contributions. This remains true for model (ii); instead, for 
model (iii), only two divergent terms remain due to a cancellation of poles which is easily 
understood by realizing that the extrinsic curvature seen from the interior and exterior has 
opposite sign. 

As the physical system we consider, as described, the classical part ( |4.7| ) and the ground 
state energy of the quantum field together, and write for the complete energy 

E = E c i ass + E vac . (4.9) 
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In this context, the renormalization can be achieved by shifting the parameters in E c i ass by 
an amount which cancels the divergent contributions and removes completely the contribution 
of the corresponding heat kernel coefficients. In detail we have 



P ~* P ~ 71 — 9 n + ln 9 > (J ^ CF + 
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(4.10) 



h + 



After the subtraction of these contributions from E vac we denote it by E™™, see eq. (pL5] 
and the complete energy becomes 



J5 



ci _|_ zjiren 
J^class i -C^uac ; 



(4.11) 



with the bare parameters in E c i ass replaced by the renormalized ones using eq. 4.10). In our 
renormalization scheme, we have defined a unique renormalized groundstate energy E^. 

Now, having discussed in detail the subtraction procedure and the structure of diver- 
gences, let us come to a full evaluation of E vac . The equations of Section |2| serve as a starting 
point. Choosing N = 3 in eq. ( |2.13[) , Z{s) is finite at s = —1/2 and can be used for the 
numerical evaluation of E vac for (in principle) arbitrary mass m. For Aj,(s), i = —1, 0, 1, 2, 3, 
a representation valid for \ma\ < 1/2 has been given, see eqs. ( |2.25|) , (|2.26|) and ( 2.2S| ), but 
here we need a representation of Ai(s) beyond that range. This is a purely technical compli- 
cation explained in Appendix |X[ The final result for A-i(s) and Aq(s) about s = —1/2 is 
given in eqs. (A. 9) and (|A.10| ), The remaining Aj(s) have the form, see eq. ( 2.2^ ), 
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,^^^E^±^ /(s;1+2c;c+i/2) , 

T(s) ^ Q {ma) t+Zc T(c + t/2) 



(4.12) 



with 



f(s;c;b)= Yl vC h + 

^=1/2,3/2,... V 




-s-b 



(4.13) 



The remaining task is to calculate /(s; c; b) for the relevant values of c and b about 
s = —1/2. This is a systematic calculation sketched also in Appendix |A|. It is slightly 
simplified by realizing the recurrence 



f(s- c; b) = (ma) 2 [f(s; c - 2; b - 1) - f(s; c - 2; b)] . 



(4.14) 



In summary, all analytic expressions needed for the numerical evaluation have been pro- 
vided. 

The result for the Casimir energy for the interior region of the ball is shown in Fig. |3| for 
a = 1 as a function of m. 

For very small values of the argument, ma, the function takes negative values, whereas for 
larger values it is positive and tends to zero for (ma) — ► do. Thus, as a function of the mass, 
attractive as well as repulsive forces are possible and the influence of the mass is by no means 
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Figure 3: Plot of the renormalized vacuum energy E. 
the radius. 
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measured in units of the inverse of 



negligeable. This is different from the example of parallel plates [21, 68], where for large mass 
the influence of the mass is exponentially damped away and is of a very short range. The 
origin of this different behaviour is intimately connected with the extrinsic curvature being 
non- vanishing or vanishing. This is already seen from the heat-kernel expansion ( f4.1[ ), where 
the leading contribution to E%™ is proportional to a 5 / 2 (— A)/m. For the ball, this coefficient 
is non- vanishing; however, for parallel plates this one as well as all the following coefficients 
are vanishing, what explains the very different behaviour of the renormalized Casimir energy 
as a function of the mass for this two configurations. 



4.2 Scalar field in the presence of an external background field 

In the above described situation, "empty space" is supplemented by boundary conditions. 
However, clearly (perfect) boundaries are an idealization and no boundary made of matter 
can be perfectly smooth. Seen as modeling some distribution of matter which interacts 
with the quantum field, it might be more justified to introduce an external potential and 
to quantize a theory within this potential. This will change the quantum modes and the 
associated spectrum and, as a result, the energy of the vacuum. For certain types of external 
confinement this situation has been called Casimir effect with soft or semihard boundaries 
|129| , |130fl . As mentioned, apart from this type of setting, in many situations the potential is 
provided by classical solutions to nonlinear field equations. 

Again, we will start the consideration by describing our concrete model and its renor- 
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malization, introducing also various notations used in the following. We will consider the 
Lagrangian 

L = -$(□ — M 2 — A$ 2 )$ + -<p(U - m 2 - A'$ 2 )v?. (4.15) 

Here, the field $ is a classical background field. By means of 

V(x) = A'$ 2 (4.16) 

it defines the potential in ( 4.15| ) for the field <^(x), which should be quantized in the back- 
ground of V(x). As explained below, the embedding into this external system is necessary in 
order to guarantee the renormalizability of the ground state energy. Actually, this is clear al- 
ready from the beginning, because the external system needs to comprise of the counterterms 
of a (A(^ 4 )-theory. 

The complete energy 

E[*]=E da8a [*] + E vac [*] (4.17) 

of the system consists, as before, of the classical part and the contributions resulting from the 
ground state energy of the quantum field ip in the background of the field <&. The classical 
part reads 

E das M = ^V g + X -M 2 V X + XV 2 , (4.18) 

with the definitions V g = J d 3 x(V$) 2 , V\ = J d 3 x$ 2 and V 2 = J cfix® 4 . Here M 2 and A are 
the bare mass, respectively coupling constant, which need renormalization. Continuing as in 
Section ^j, for the ground state energy one defines |28j 

EvaM = \ £(A 2 n) + m 2 ) 1 /^, (4.19) 

(n) 

where \x is the arbitrary mass parameter and s is the regularization parameter, which has to be 
put to zero after renormalization. Furthermore, A( n ) are the eigenvalues of the corresponding 
Laplace equation 

(_A + V(x))4> (n) (x) = \ 2 {n) 4> in) (x). (4.20) 

For the moment, we assume the space to be a large ball of radius R as an intermediate 
step to have discrete eigenvalues and, thus, a discrete multiindex (n) and, in addition, to 
avoid pure volume divergences. For simplicity we impose Dirichlet boundary conditions and 
explain below under which conditions on the potential V(x) the final result does not receive 
boundary contributions. 

Expressing the ground state energy fl4.19|) in terms of the zeta function 

Cy(s)=£(A 2 n) + m 2 r (4-21) 

(n) 

of the wave operator with potential V(x) as defined in ( |4.20| ), one has 

EvaM = 1 -Cv{s-1/2)li 2s . (4.22) 
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The residue of Cv(~ 1/2) is determined by the ai heat-kernel coefficient associated with 



the spectrum A 2 n ^ + m 2 , and the classical energy ([4.18 ) provides another example for our 



previous viewpoint, that it should consist of the terms contained in a%. In order to fix the 
finite renormalizations, we impose again the normalization condition 

Jim^pHO, (4.23) 
which is implemented along the lines leading to eq. (4.6). For convenience we state it again, 
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(4.24) 
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2 02, 



with aj, i 



0, 1/2, 2 the heat-kernel coefficients of the Laplace equation ( 4.2C| ). This is a 
good moment to consider the limiting behaviour for R — > cxd of the boundary terms. Before 
actually performing the limit R — > cxd we must clearly state what the ground state energy 
in the given context is meant to be. It is a quantity that describes the influence of the 
background potential V(x) on the vacuum energy and, as such, is to be compared with the 
field free case V(x) = 0. Thus, from definition ( 4.19Q , we will subtract the Casimir energy 
of a free field inside a large ball of radius R in order to normalize E vac [Q = 0] = 0. In 
eq. (4.24), this subtraction cancels the contributions of ao and ai/2> i n a i it cancels the 
potential independent boundary terms. In 03/2 the term proportional to V(R)R 2 survives, 
see eq. (|3.10| ), and, in the limit R — > 00 vanishes only if V(r) ~ r~ 2 ~ € , e > 0, for r — ► 00. As 
is easily seen by dimensional arguments, the boundary contributions of the higher coefficients 
depending on V(r) will then vanish too. 

In summary, under the assumption that V(r) ~ r~ 2 ~ e for r — > 00, in the limit R — » cxd 



no boundary contributions will appear and the normalization condition ( 4.23 ) is achieved by 
a renormalization of the mass M of the background field, 



M — 3 

and the coupling constant A by 

A - 
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(4.25) 



(4.26) 



This is an immediate result of the volume contributions to the heat-kernel coefficients, 
ai = -(4tt)~ 3 / 2 / ci 3 xy(x) and a 2 = (l/2)(4vr)- 3 / 2 / d 3 xV 2 (x). The kinetic term F 9 in E[$] 
suffers no renormalization. 

By defining 



div 



(4.27) 



one obtains the finite groundstate energy, which is normalized in a way that the functional 
dependence on <I> 2 present in the classical energy is now absent in the quantum corrections 



-'--'vac ' 
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Let us note that this is just the well known general renormalization scheme written down 
here explicitly in the notations needed in our case. 

For the calculation of i£ mc [<&] let us take further advantage of the background field 3>(r) 
being spherically symmetric. The multiindex (n) — * n,l,m consists of the main quantum 
number n, the angular momentum number I and the magnetic quantum number m. In polar 
coordinates the ansatz for a solution of the wave equation (4.20) reads 



1 

= -4>n,l(r)Y lm {6,ip) 



where the radial wave equation takes the form 



dr 2 



1(1 + 1) 



+ A 2 



n. i 



<f>n,l(r) = 0. 



(4.28) 



(4.29) 



Now we use the standard scattering theory within r £ [0, oo) and take the momentum p 
instead of the discrete A n) j. Let 4> p ,i(r) be the so called regular solution which is defined so 



as to have the same behaviour at r — > as the solution without potential 



t>p,l\ 



jlipr) 



(4.30) 



r-»0 



with ji the spherical Bessel function [131]. This regular solution defines the Jost function fi 
through its asymptotics as r — ► oo, 



<t>i, P { r ) 



fl(p)h, (pr) - fi(p)hf(pr) 



(4.31) 



where h^(pr) and h[(pr) are the Riccati-Hankel functions [131]. The analytic properties of 
the Jost function flip) strongly depend on the properties of the potential V(r). If in addition 
to V(r) ~ r ~ 2 - e fo r r — > oo ; we impose V(r) ~ r _2 + e for e — > and continuity of V(r) in 
< r < oo (except maybe at a finite number of finite discontinuities), one may show that 
the Jost function is an analytic function of p for Qp > 0. It vanishes in the finite set of points 
p = iK n i of the positive imaginary half axis, corresponding to the bound states with energy 

2 

K n,l- 



Now we use the Jost function to transform the frequency sum in eq. (4.21) in a contour 
integral. Let us assume that the support of the potential is contained in the cavity of radius 
R. Then the above eq. ( 4.31| ) gets exact at r = R and may be interpreted as an implicit 
equation for the eigenvalues p = \ n j. Choosing Dirichlet boundary conditions at r = R, 
(j)p t l(R) = 0, it reads explicitly, 



fi(p)hJ-(pR)-f?( P )h+(jpR) = 0. 



(4.32) 



As already mentioned, ultimately we are interested in the limit R — > oo and in that limit 
the results will not receive boundary contributions, once we assume that V(r) ~ r ^ 2 ~ e fo r 
r — > oo. 

Let us now consider the ground state energy associated with the eigenvalues determined 
by (4.32). As described in detail in Section ||, we represent the frequency sum in ( [4.19[ ) by a 
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contour integral. Using eq. 



), one immediately finds 
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+ M 2s £a + V2)£(™ 2 -<*) 1/2 - 

1=0 n 



(4.33) 



with 



the energy eigenvalues of the bound states with given orbital momentum I. The 



denominator in the logarithm provides the subtraction of the Minkowski space contribution, 
where fi(p) = 1. The contour 7 is chosen counterclockwise enclosing all real solutions of 
eq. ( |4.32p on the positive real axis. In the limit of the infinite space the negative eigenvalues 
become the usual boundstates and the \ n j > turn into the scattering states. 

For the calculation of ( [4.33 ), as the next step, one deforms the contour 7 to the imaginary 
axis. A contour coming from ioo + e, crossing the imaginary axis at some positive value smaller 
than the smallest K n and going to ioo — e results first. Shifting the contour over the bound 
state values K n , which are the zeroes of the Jost function on the imaginary axis, the bound 
state contributions in eq. (4.33) are cancelled and in the limit R — > 00 one finds 
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(4.34) 



For details of the contour deformation see Fig. ||. 

In the first step of the deformation one makes use of the following properties [131 1 , 



M-p) 
hH-z] 



-l) l hf{z). 



This is the representation of the ground state energy (and by means of ( [4.22j ) as well of 
the zeta function) in terms of the Jost function, which is the starting point of our following 
analysis. It has the nice property, that the dependence on the bound states is not present 
explicitly, being however contained in the Jost function by its properties on the positive 
imaginary axis. It is connected with the more conventional representations by means of the 
analytic properties of the Jost function. Expressing them by the dispersion relation 
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where 5i(q) is the scattering phase, we obtain from (4.34) 




E 



m 



l-2s 



-l-2s^ 



777/ 



1 - 2s 



IT 



T l-2; 



5i{q) 



(4.35) 



which gives the expression of the ground state energy through the scattering phase. From 
here one can pass to the representation through the mode density by integrating by parts. 
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Figure 4: Deformation of the contour 7 



Note that this representation can be obtained also directly from ( [4.33| ) in the limit R —* 00 
by deforming the contour 7. 



Let us add a discussion on the sign of the ground state energy. In ( |4.35[) the first contri- 
bution results from the bound states and is completely negative. The second contribution, 
which contains the scattering phase 5i(q), is positive (negative) for an attractive (repulsive) 
potential, i.e., for V(r) < (V(r) > 0) for all r |131[ |. So the regularized (still not renormal- 
ized) ground state energy E vac [$] ( f4.35| ) is positive for a potential which is repulsive for all r 
(there are no bound states in this case) and it is negative for a potential which is attractive 
for all r. Now, if we perform the renormalization in accordance with ( 4.25| ) and ( f4,26| ), we 
obtain 

rpren p 

^vac ~ ^vac 
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+— [ - +log-^ - 2 1 / drr z V(r) z . (4.36) 







This expression is finite for s — > 0, i.e., when removing the regularization. But due to 
the subtracted terms there is no longer any definite result on the sign. Note that this is in 
contrast to the case of a one-dimensional potential, where it had been possible to express the 
subtracted terms through the scattering phase ]56|| . 

Let us continue with a detailed analysis of the ground state energy, eq. ( [4.34 ). As we 



have easily seen using heat-kernel techniques, the non-renormalized vacuum energy £J mc [$] 
contains divergencies in s = 0, see eq. ( (4.24| ), which are removed by the renormalization 
prescription given in eqs. ( |4,25 ) - fl4,27| ), resp. Q4.36] ), The poles present are by no means 



obvious in the representation ( 4.34| ) of E vac [$>]. However, in order to actually perform the 



renormalization, eq. ( 4.27|) , it is necessary to represent the groundstate energy eq. ( 4.34j ) in 



a form which makes the explicit subtraction of the divergencies visible. This will be our first 
task. 

As is known from general zeta function theory as well as one sees from simply counting 
the large momentum behaviour, the representation eq. ( [4.34| ) of i? m c[^] will be convergent 
for Res > 2. However, for the calculation of the ground state energy we need the value of 
eq. ( 4.34j ) at s = 0; thus, an analytical continuation to the left has to be constructed. The 



basic idea is the same as before: adding and subtracting the leading uniform asymptotics of 
the integrand in eq. (fO§) f72|. Let 

E vac [<f>]=E f + E as , (4.37) 

where 

E f = - C -^M^ £)(J + 1/2) jdk [k 2 - m^-^ilnhiik) - Intf"*^*)] (4.38) 



and 
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As described in detail, the idea is that as many asymptotic terms are subtracted as 
necessary to take s = in the integrand of Ef. This term will then (in general) be evaluated 
numerically. 

In E as the analytic continuation to s = can be done explicitly, showing that the pole 
contributions cancel when subtracting E^ c [^], eq. ( [4.24 ). 



The first task, thus, is to obtain the asymptotics of the Jost functions. Contrary to 
the case with boundary conditions, the asymptotics cannot be taken just from tables, but 
need to be calculated itself. Fortunately, this may be done by using the integral equation 



(Lippmann-Schwinger equation) known from scattering theory [131]. For the Jost function 
one has (y = I + 1/2) 



00 



fl(ik) = 1+ / dr rV{r)4>, ik {r)K v {kr), (4.40) 
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with the regular solution given by the integral equation 

r 

faM r ) = /,(A:r)+|dr / rM/,(A ; r)^(A;r , )-/,(A;r / )^(fcr)]F(r')0 Mfe (r / ). (4.41) 
Q 

General zeta function theory tells us that the divergence at s = contains at most terms 
of order V 2 . Thus, one might expand ha. fi(ik) in powers of V and take into account only the 
asymptotics of terms up to 0(V 2 ). The expansion in powers of V is easily obtained. Using 
eqs. ( [4.40D and ( |4.41[) one finds 



dr rV(r)K u (kr)I u (kr) 
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(4.42) 



Now, the uniform asymptotics for I — > oo of In fi(ik) is essentially reduced to the well 
known uniform asymptotics of the modified Bessel functions K v and I u , (9.7.7) and (9.7.8) 
in ||J. With the notation t = + {kr/v) 2 and rj(k) = y/1 + (kr/f) 2 + ln[(kr/u)/(l + 

y/1 + (kr/u)' 2 )], one finds for v — ► oo, k — » oo with kjv fixed, 
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Using these terms in the rhs of eq. ( [4.42| ) we define 
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(4.43) 



Thereby the r'-integration in the term quadratic in V has been performed by the saddle- 
point method using the monotony of f](k). Now, by means of ( 4.43| ) the limit s — > can be 
performed in eq. ( [4. 38 ) and we obtain 

i 00 °f , fi 

E f = -- £(i + 1/2) y ^ Vfe2- m 2— (ln/,(tfc) - ln/f S2/m (^)) , (4.44) 

Z = 777, 
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a form which is suited for a numerical evaluation. 

For E as at s = one might explicitly find the analytical continuation. First of all, the 
/c-integrals may be done using 
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ok 
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This naturally leads to the functions encountered already in the case of boundary condi- 
tions, 



f(s;c;b)= J2 vC 1+ — 

i/=l/2,3/2,... V 
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(4.45) 



In terms of these functions, for E as we obtain 
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The relevant expansions about s = of the f(s;c;b) are found in Appendix |A]. All 
divergences are made explicit, and 



takes the compact form Jf. 
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(4.46) 
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valid for any potential with the above mentioned properties and very suitable for numerical 
evaluation. 

The remaining task for the analysis of the renormalized ground state energy 
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in the presence of a spherically symmetric potential is the numerical analysis of Et. To 
achieve that, a (as a rule numerical) knowledge of the Jost function fi(ik) is necessary. This 
will be obtained from a numerical knowledge of the regular solution at a single point [73|. 

Under the assumption that the potential has compact support, let's say V(r) = for 
r > R, the regular solution may be written as 



<M0 = u l>p (r)e(R - r) + - fi(p)hf(pr) - /,* (p)h+(pr) G(r - R) 



(4.47) 



Assuming continuity of the regular solution and its derivative, the matching conditions 



arc 



u lA R ) 
u l, P ( R ) 



fi(p)K(pR)-fi(p)K(p R \ > 
~Mp)hy'( P R) - ft{p)hf'{pR) . 

These conditions provide already the Jost function in terms of the regular solutions, 

flip) = ~ (put, P (R)ht'(pR) - u' hp {R)hf(jpR)) , (4.48) 



where we used that the Wronskian determinant of hi is 2i. 



Compared to integral representations of the Jost function [ 131 



fl(ik) = 1 + / dr rV(r)^ lk (r)K u (kr), 
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(4.49) 



eq. ( 4.48| ) has the advantage that the solution is only needed at the point r = R. This 
simplifies the numerical procedure considerably. 

In order to determine a unique solution of the differential equation ( [4.29 ), we need to 
pose an initial value problem. Given that (fii tP (r) is the regular solution, we have for r — > 



,(r) = ui >p (r) ~ ji(pr) 
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The natural ansatz thus is 
ui, P (r) 



pr 
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T(/ + 3/2) V 2 

with the inital value gi :P (0) = 1. The differential equation for gi tP (r) reads 

l^ + 2—^-V(r)+pAg Wp (r) = 0, 
I ar z r dr I 

or, going to the needed imaginary p-axis and writing it as a first order differential equation 
with (d/dr)g^ ip (r) = u/,j p (r), 



/ 9i,i P (r) 
dr I vi tip (r) 



1 W 9i,ip{r) 
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(4.50) 
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To fix the solution uniquely we only need to fix vn p (0). A power series ansatz for gn p (r) 
about r = shows that, for V(r) = 0(r~ 1+e ), the condition reads v/ j j p (0) = 0. With this 



unique solution of (4.50), the Jost function takes the form 

2 fvR\ l+3 / 2 ( 11 
fl{iT>) = r(/ + 3/2) \T) {9i,ip( R ) K l+s/2(.pR) + -gi P (R)K l+1/2 (pR)j . (4.51) 

Finally, doing a partial integration and the substitution q = \Jk 2 — m 2 in eq. (fy||, the 
starting point for the numerical evaluation used for Ef is 
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Now we are prepared to use eqs. ( [4.46| ) and ( 4.52 ) for the calculation of E r v ™. We choose 



a potential with compact support <3?(r > R) = 0. For the numerical analysis to come, we 
have used the dimensionless quantities 



r 



e = E vac R, p = mR, P = ^ 
V(r) = X^ 2 (r) = ^(p). 



The example we consider here is |7S 

16ap 2 (l - p) 2 



<P(P) 



a + (l-2pf 



which is a kind of a spherical wall; the parameter a allows to vary the shape of the potential 
(see Fig. |5|). 

The most determining property of the potential V(r) is the number and depth of bound- 
states. The parameter dependence on A' of the / = boundstates, determined as zeroes of 
the Jost function, are given in Fig. || As is clear, number as well as depth of the boundstates 
increases with increasing (—A'), and, as is seen, with increasing a. 

To show clearly the way this influences the vacuum energy, we consider the dependence 
of the energy for fixed A' as a function of the parameter a, Fig. ^, as well as for fixed a = 1 
as a function of A', Fig. 0. 

Let us start with a description of Fig. ^. For a large enough, the contribution of the 
bound state(s) to the vacuum energy is large enough to overcompensate the scattering state 
contributions and the energy is negative. At some "critical value" of a, the energy of the only 
remaining bound state is so small, that the positive contributions of the scattering states get 
the upper hand and the vacuum energy is positive. Finally, for a — > we have normalized 
Eyaci^ = 0] = 0) such that at some point, starts to decrease again with decreasing a. 

The same features are clearly recovered in Fig. |?]. For large enough (—A') a negative 
vacuum energy is obtained. With decreasing (—A') the energy increases and can be positive 
or negative depending on the parameter p; for fixed value of R this means the mass to. Again, 
at some point the scattering states get the upper hand, the energy gets positive and tends to 
zero as A' tends to zero. For positive A' the energy starts to increase again, being virtually 
identical for A = ±2, and it seems that -f^ac i s a monotonically increasing function of A' for 
A' > 0. 
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Figure 5: Energy of bound states for several shapes of the potential and negative A'. (The 
inset shows tp(p).) 



a= 10 (2 BS) 
a= 0.1 (1 BS) 
a= 0.001 (1 BS) 
a=0.0001 (no BS) 
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Figure 7: Vacuum energy for different magnitudes A' with equal shape parameter a = 1. 



Exactly the same features are observed if instead of the chosen spherical wall a lump- 
like potential concentrated around r = is chosen, and it is expected that these features 
hold for any kind of potential |7ql . In addition, the numerical analysis seems to indicate 
that, without bound states, the vacuum energy is always positive, indicating that scattering 
states contribute positively to the vacuum energy (a property we actually used already in the 
description of the figures), a result for which no analytical proof exists in three dimensions. 



5 Conclusions 

In this contribution we have provided methods for the analysis of zeta functions associated 
with second order differential operators, when spherically symmetric external conditions are 
present. Basic analytical skills applied are contour and Merlin-Barnes integral representations 
and, at a second stage, the Barnes zeta function. Due to the close connection between 
zeta functions and various other spectral functions as the heat-kernel and determinants, 
applications to the calculation of heat-kernel coefficients and ground state energies have been 
given. 

The determination of heat-kernel coefficients is achieved by a combination of different 
methods, see Section ||. As a first step write down the general form of the coefficients in terms 
of geometrical invariants. The possible invariants strongly depend on the boundary conditions 
considered as has been briefly described || 123]. Next one can consider the special case 



calculations of the ball given in Section 3.2 with this general form and fit unknown numerical 



constants, as has been explained in detail in Section 3.3. Together with information obtained 



from the product formula, this gives a very good starting point to apply the functorial 



techniques developed in []96[. Experience shows that, seemingly, any boundary condition 



can be successfully attacked in that way. 
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Afterwards, we concentrated on the analysis of vacuum energies. For the case of a massive 
scalar field inside a spherical shell we calculated the Casimir energy for the case when the field 
satisfies Dirichlet boundary conditions. Contrary to what one expects from the example of 
parallel plates, the energy depends strongly on the mass and, as a function of it, attraction as 
well as repulsion is possible. Mathematically, the calculation of the vacuum energy amounts 
to an evaluation of the associated zeta function at the specific point s = —1/2. With the 
help of the analytical approach provided in Section ^ this reduces to a numerical analysis of 
convergent integrals. 

If an external potential is present we have expressed E^ through the associated quantum 
mechanical scattering dates, namely the Jost function. This is explicitly calculated by solving 
numerically a first order linear differential equation and by performing numerical integrations. 
The vacuum energy behaves in a way one can understand from physical reasons and the 
definition introduced by the normalization condition (|4.23| ) seems to be very reasonable. 

Let us emphasize that the external potential can be a solution of some classical nonlinear 
equation. But to apply the method one needs just any spherically symmetric potential of 
otherwise whatever shape to come immediately to a determination of E^. This is hoped to 
find further applications in the context of the fields mentioned briefly in the Introduction. 
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A Representations for the asymptotic contributions 



In this Appendix we derive explicit representations of the asymptotic contributions in the 
Casimir and ground state energies for massive fields. Let us start with A-\(s) for the massive 
scalar field, eq. ( 2.15| ), which, a little bit more explicit reads 
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of which we need the analytical continuation to s = —1/2. With the substitution t 
ixujd) 2 — m 2 , this expression results into the following one 
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where the Mellin integral representation for the single factors has been used. As we see, 
the /3-integral is well defined. Introducing a regularization parameter 5, A-i(s) can then be 
written as 



A_ 1 ( S ) = lim Al l (s,S) + Al 1 (s,5) 

0— >U L 
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Let us proceed with the remaining pieces. In A_^(s, S) two of the integrals can be done, 
yielding 
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For J 4?_ 1 (s,(5), one gets 
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Adding up (|A.4j) and (|A.5|) yields 
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where in the last step a partial integration has been performed such that the x-integral is 
well behaved at s = —1/2. 

This is a form suited for the treatment of the angular momentum sum, which is performed 
using 
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where e — > is understood and appropriate analytic properties of the function f(u) are 
assumed. This is often used in finite temperature field theory and is equivalent there to 
switching from imaginary time to real time [132 1 . When expanding the function f(y) in a 
Taylor series, one arrives at the well known Euler-Maclaurin summation formula (a thorough 
treatment of the Euler-Maclaurin summation formula can be found in Ref. 133]). 
For A_\{s) the relevant application of eq. (|A.7|) is 
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the following final form of A_i(s) is obtained, 
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In the same way, for 
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it follows 
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For the remaining asymptotic contributions, the central spectral function is, see eq. (4.13), 



f(s;c;b) = »° ) 



-s-b 



To deal with all values of c and b needed, in addition to eq. flA.Sj ), we need 

2 V S ir(n+l/2)r( S -n-l/2) ^ 2n+1 



E » 2n 1 + 

v=l/2,3/2,... \ 



2 



(A.ll) 



Using partial integrations one can obtain representations valid for values of s needed for 
the asymptotic contributions. To simplify notation, in the following we shall use x for (ma) 
and /(c; b) = /(— 1/2; c; b). A full list of ingredients needed is the following: 

1/2) = -I^ + i, 



d_ 

ds 



_ 1/2 /( S ;l;l/2) = --x 2 -2 



dv 



1 + e 2 ™ 
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/(l; 3/2) = — — — + x 2 lnx + x 2 / 
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/(1;1) =2^ / ^ 



/ d 1 



Vdi/ 1 + e 2 ™ 
/ d 1 



/(l;2) = -2^f rf ,(- I 



/(3;2) = -x 4 + 2x 4 / <fr 
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+ e 



did v 2 



d 1 
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1I/2 
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+ e 
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x dv 
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1 
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df I 1/ df ^ dv 1 + e 27rl/ 
+ (lnx- ll/12)x 8 

d ( I d I d I d v e 

dv 
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2 2 
^ — x 



dv \ v dv v dv v dv 1 + e 27W 
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2 2 
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